Anomalous low-field classical magnetoresistance in two dimensions 
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The magnetoresistance of classical two-dimensional elec- 
trons scattered by randomly distributed impurities is investi- 
gated by numerical simulation. At low magnetic fields we find 
for the first time a negative magnetoresistance proportional 
to \B\. This unexpected behavior is shown to be due to a 
memory effect specific for backscattering events, which was 
not considered previously. 

PACS numbers: 05.60. +w, 73.40.-c, 73.50.Jt 



The problem of magnetoresistance in metals and semi- 
conductors has a long history, and a large amount of 
work, both experimental and theoretical, was devoted 
to this subject. Theoretically, the simplest situation is 
that of two-dimensional degenerate non-interacting elec- 
trons moving in a plane perpendicular to the magnetic 
field and scattered by static impurity potential (the 2D 
electron gas at low temperatures is also the situation for 
which most of the experiments were done). In this case 
the relevant electron energy is equal to the Fermi energy, 
and the conventional Boltzmann-Drude approach yields 
zero magnetoresistance, i.e. the longitudinal resistivity, 
Pxx^ does not depend on magnetic field. Thus the origin 
of the observed magnetoresistance should be looked for 
beyond the Boltzmann theory. 

Most of the research in this domain was focused on the 
low-field magnetoresistance arising from quantum inter- 
ference effects (weak localization) ^ . Although the pio- 
neering work of of Baskin et al |^ long ago demonstrated 
the importance of non-Boltzmann classical memory ef- 
fects in magnetotransport (see also Refs. [^,||), only in 
recent years it was fully recognized that there are not 
only quantum, but also purely classical reasons for mag- 
netoresistance, 11,^,0,^1), which may be either negative 
or positive depending on the type of impurity potential. 
Whatever is the case, the classical magnetoresistance ap- 
pears as a consequence of memory effects which are be- 
yond the Boltzmann approach. 

For the case of strong short-range scatterers, with 
which we are concerned, there is a large negative mag- 
netoresistance at /3 > 1, (Hill] where (3 = WcT, lOc is 
the electron cyclotron frequency, and r is the momen- 
tum relaxation time. This is related to the existence 
of "circling" electrons, |^,^ which never collide with the 
short-range scattering centers, the fraction of such elec- 



trons being 

P exp(-27ri?/^) exp(-27r//3), (1) 

where R = v /uoc is the cyclotron radius, v is the electron 
(Fermi) velocity, and £ — vt ~ {Nd)~^ is the electron 
mean free path. For the case of isotropic scattering, the 
corresponding magnetoresistance is given by the simple 
formula (valid with an accuracy better than 2%): 

Pxx = Pq{1-P), (2) 

where po is the zero-field resistivity. This formula is in 
good agreement with numerical simulation |9j in the limit 
d/£ = Ncf — > 0, where d is the effective diameter of 
the scattering centers and N is their concentration. Eq. 
(2), which is valid in this limit, predicts an exponentially 
small magnetoresistance at low fields, when (3 << 1. 
However for finite values of the parameter d/£ a parabolic 
negative magnetoresistance was found at /3 << 1. |^ 

The purpose of this Letter is to report a numerical 
study of the low-field classical magnetoresistance of a 2D 
electron gas with short-range scattering centers. Unex- 
pectedly, our numerical simulation reveals a new charac- 
teristic magnetic field defined by the relation (3 = d/ £ « 
1. For /3 << d/£ the (negative) magnetoresistance is lin- 
ear in B, while outside this interval, for d/£ « (3 « 1 
the dependence is parabolic. We explain the physics of 
this previously unknown phenomenon and we discuss its 
relevance to the experimental results for magnetotrans- 
port in random antidot arrays. |[lO| , pT[ | 

In our simulation a point particle (electron) with a 
given absolute value of its velocity, v, is scattered by 
disks of diameter d randomly positioned on a plane inside 
a square box with periodic boundary conditions (the box 
size is a thousand disk diameters). Both the hard-disk 
(Lorentz) model, which exhibits anisotropic scattering, 
and a modified model with isotropic scattering are stud- 
ied. To characterize the coverage, we introduce a dimen- 
sionless concentration of scatterers c = iiNd^ /A << 1. 

The simulation procedure is identical with the one de- 
scribed in Ref. |^ : an initial electron position is chosen at 
random with an initial velocity along the x direction, the 
successive velocity directions are determined after each 
collision and the conductivity tensor is determined by 
calculating the integral of the velocity-velocity correla- 
tion function over a time oi t = 50t. For each value of 
field and disk concentration we take an average over 10^ 
independent disk configurations and, for each configura- 
tion, over lO'' independent trials for the initial electron 
position (10 time more than in ref. Q). 
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Our main result is presented in Fig. 1, which shows 
a characteristic anomaly in magnetoresistance, Ap = 
Pxx{B) — pq, around zero magnetic field followed by 
the anticipated parabolic dependence on f3 (hard disk 
model) . The anomaly is more visible when Ap is plotted 
against /3^, see inset. For still higher fields the curve fol- 
lows Eq. 2, as described in Ref. ||^. Similar results were 
obtained for isotropic scattering and other values of c. 
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FIG. 1. Numerical results for tlie magnetoresistance Ap/po 
as a function of /3 = uJcT for c = 0.1 (hard disk model). Inset: 
same data plotted as a function of /3^. Here and below the 
error bars are within the symbol size. 

To understand this result we need to discuss in some 
detail the nature of the classical non-Boltzmann correc- 
tions to the Drude resistivity. The Boltzmann approach 
is equivalent to a random redistribution of all scatterers 
after each collision. In reality there are memory effects of 
two types: (i) the particle may re-collide with the same 
impurity and (ii) the particle's trajectory may repeatedly 
pass through a space region which is free of impurities 
(the effect of circling electrons is of this type). 

Consider first the case B = 0. The classical corrections 
to the Boltzmann equation due to re-collisions with the 
same scatterer (i) were studied long ago. These pro- 
cesses are responsible for the 1 /i^ tail in the velocity cor- 
relation function ||l^ and for the increase of the resistiv- 
ity compared to its Drude value. The relative correction 
for small c is proportional to cln(l/c) (or equivalently, to 
{d/£) \n{£/d)), the logarithmic term being due to the sim- 
plest re-collision process 1 — > 2 ^ 1, while longer return 
loops give a smaller correction on the order of c. 

As to the type (ii) processes, their role for B = was 
not, to our knowledge, well understood so far. Consider 
a particle which, after going a distance x » d without 
collisions, experiences backscattering at an angle </> = tt 
and then returns to the initial point (Fig. 2a). The 
probability of this round trip of length 2x is proportional 
to exp{—x/£), not to exp{—2x/£), as would suggest the 
conventional approach, since the existence of a free cor- 
ridor of width d allowing the first part of the journey 
guarantees a collision-less return. This is not the case for 
scattering angles outside the interval on the order oi d/x 
around the value <j} = t:, when the probabilities for a free 



path X before and after collision become independent and 
equal to exp(— x/^) each (Fig. 26). 

Since typically x ^ £, the probability of backscattering 
in the interval Ac/) ^ d/£ around (f) — n is enhanced 
(the reversed velocity is conserved for a longer time), and 
this should lead to an additional increase of resistivity 
on the order oi d/£ ~ c, i.e. same order of magnitude 
as the contribution of return loops involving two or more 
intermediate scattering. One can say that the existence 
of a free corridor effectively enhances backscattering in 
the interval A0, roughly by a factor of 2, thus increasing 
the transport cross-section by an amount ~ dA(j). 

We attribute the low field anomaly in Fig. 1 to the 
influence of magnetic field on this effect. In the presence 
of even a small magnetic field the electron trajectories 
before and after collisions can not any more follow the 
same path (Fig. 2c). At high enough fields this kills 
the memory effect and, as a consequence, reduces the 
resistivity (Fig. 2d). Thus a negative magnetoresistance 
appears with a characteristic magnetic field defined by 
the relation (3 = d/£ << 1. 
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FIG. 2. Illustration of the memory effect in backscattering. 
a - Backscattering at angle vr, B = 0. The particle returns 
through the same corridor, b - Scattering at an angle far from 
TT, B = 0. The two corridors are different and the probabilities 
of free flights before and after collision are uncorrelated. c - 
Scattering at angle tt in low magnetic field, R >> /d. The 
overlap of the two corridors diminishes, d - Same as c, but 
for higher field, R « /d. The corridors cease to overlap. 

We do not have a regular method for calculating an- 
alytically the magnetoresistance at low fields. However 
some qualitative conclusions may be drawn from simple 
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geometric considerations presented in Fig. 2. 

Consider again backscattering by an angle exactly 
equal to tt but in the presence of magnetic field (Fig. 
2c) . In order to have collision-less paths of length x be- 
fore and after scattering, the centers of all disks should be 
outside the corridor of width d surrounding these paths. 
The probability of this is proportional to P = exp{—NS), 
where N is the disk concentration and S is the joint area 
of the two corridors. The overlapping region should not 
be counted twice. 

While at B = there is full overlap (Fig. 2a), S — xd 
and P = exp(— a:/£), in the presence of magnetic field the 
overlap diminishes and the relevant area increases. In the 
low- field limit, this increase can be easily calculated to be 
AS* = x^ /{iR) ~ B, so that P decreases linearly with B. 
This means that the (negative) magnetoresistance is lin- 
ear in B for R » f' /d, or (3 « d/£. For higher fields, 
such that /3 >> d/£ the two corridors practically cease to 
overlap (Fig. 2d) and one has S — 2xd, P — exp{~2x / £) . 
Similar considerations apply to backscattering in the in- 
terval A(f> ~ d/i. 

A similar contribution comes from the influence of 
magnetic field on the probability of the simplest re- 
collision process 1 — > 2 ^ 1, which necessarily involves 
backscattering in the same angular interval A^. At 
B = the memory effect increases the relative contri- 
bution of this process to the resistivity by an amount 
on the order of c, and again the curving of the trajecto- 
ries in magnetic field will increase the total area S and 
thus reduce the probability of this process. In the low- 
field limit one finds again that the area increase, AS", 
is linear in B. These qualitative considerations lead us 
to the following conclusions concerning low-field classical 
magnetoresistance in two dimensions in the presence of 
appreciable backscattering. 

1) A new characteristic magnetic field exists, at which 
the classical parameter f] = lUcT is small: Pc — d/£ ~ 
c « 1. 

2) The total drop of resistivity in the region /3 < /3c is 
on the order of d/£ ~ c. 

3) At /3 << /3c the resistivity pxx is linear in magnetic 
field exhibiting the \B\ cusp observed in our simulation 

@- 

4) For /3c << /3 << 1 only quadratic in B correc- 
tions remain, which are on the order of c/3^. (It can be 
shown that these corrections come from the influence of 
the magnetic field on the contribution of return loops). 

This means that at low fields, /3 << 1, the magnetore- 
sistance, Ap, is described by the formula: 

Ap/po = -c{f{(3/c) + A(3'), (3) 

where A is a numerical constant and /(^) is a function 
which behaves as |^| for small values of its argument and 
saturates at some value on the order of 1 for |^| >> 1. 
A complete theory should give the explicit form of the 
function /(^) and the numerical value of the constant A. 
For the hard disk model, we find /(^) « 0.04|^| for ^ < 2, 
and A « 0.3. 



Note, that Eq. 3 implies that at low field the slope 
of Ap/po as a function of f3 is independent of c and that 
the derivative {d/df3){Ap/ po) has a jump at B = which 
is a numerical constant depending only on the relative 
efficiency of backscattering (« 0.08 for hard disks). 

Since the conclusions above are not derived rigorously, 
we test them by plotting the relative magnetoresistance 
in units of c as a function of the dimensionless magnetic 
field /3/c for several values of c. As follows from Eq. 3, in 
these units the magnetoresistance curves should coincide 
at low fields (where the function /(^) is linear), but di- 
verge for /3 > c. This prediction is in excellent agreement 
with the simulation results, as presented in Fig. 3, which 
convinces us that our understanding is basically correct. 
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FIG. 3. Simulation results for several values of c plotted 
in units Ap/(cpo) versus /3/c showing universal behavior at 
low fields as predicted by Eq. 3. Open circles, filled circles 
and open triangles correspond to c = 0.05, 0.1, and 0.15, 
respectively. 

Finally, we compare our simulation with the experi- 
mental results for magnetoresistance of 2D electrons in 
a disordered array of antidots. [ p^[pl| ] In these experi- 
ments the antidot diameter, d, is much greater than the 
Fermi wavelength, so that the classical approach is jus- 
tified. Our numerical simulation exactly corresponds to 
the experimental situation at low temperatures. 

Negative linear magnetoresistance at low fields, first 
observed in Ref. [p^ , quantitatively agrees with our re- 
sults as seen in Fig. 4. The necessary data being given 
in Ref. ||l^, this comparison is done without any fit- 
ting parameters. (The main uncertainty is in the ef- 
fective antidot diameter which is given in Ref. [|lO| as 
d — 0.2 — 0.3pm. We take d = 0.25pm. In order to 
re-plot the experimental curve we calculate c « 0.1 and 
/3 = 12.5B(T)). It should be noted that while there is 
an excellent agreement for (3 < 1, strong deviations are 
found at /3 > 1. The discrepancy may be due to the fact 
that the antidot array used in these experiments was, in 
fact, a strongly distorted regular lattice with a substan- 
tial long-range order. 
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FIG. 4. Simulation results for magnetoresistance at c = 0.1 
(filled circles) compared to the experimental curve of Ref.[10] 
(thick line). 

In Ref. the observed anomaly is attributed to the 
quantum (weak localization) effect. Direct comparison of 
our results with the experiment is difficult because of the 
way in which the data is presented, however we find, at 
least, a qualitative agreement both in the magnitude of 
the low-field anomaly and in the value of the characteris- 
tic magnetic field at T=1.2 K. This may be an indication 
that the experimental results are (at least partly) a mani- 
festation of the classical phenomenon presented above. A 
strong temperature dependence of the low-field anomaly 
was observed in Ref. which seems to support its 

interpretation as a result of weak localization. However, 
the observed shrinking of the low-field anomaly, as well as 
of the value po, with increasing temperature could be also 
due to the reduction of the memory effect by electron- 
electron and phonon scattering. The temperature depen- 
dence of the anomalous magnetoresistance needs further 
studies. 

The artificial random antidot array is an interesting 
object, and a systematic experimental study of magne- 
toresistance in this system is highly desirable to separate 
the classical and quantum effects in magnetotransport. 

A regular analytical calculation of the classical mag- 
netoresistance at low fields remains a challenging task. 
It would also be important to understand the transition 
from classical to quantum behavior as the ratio of the 
Fermi wavelength to the scattering diameter increases. 

Interestingly, a mechanism similar to the one consid- 
ered here is responsible for the "opposition effect" |l^] 
- the increase of the brightness of the Moon and other 
planets in a small angular interval around the position 
when the Sun, the planet, and the observer lie on the 
same straight line. 

In summary, by numerical simulation we have found 
an anomalous low-field classical magnetoresistance of 2D 
electrons in the presence of strong short-range scattering. 
A new characteristic magnetic field is found, at which the 
classical parameter /3 = oJcT is small: f3c = d/£ << 1. For 
f3 << (3c the negative magnetoresistance is proportional 
to \B\, and the slope of Ap/pQ versus /3 is independent 
of the impurity concentration. We have shown that this 
phenomenon is due to a specific memory effect associated 



with backscattering. 
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